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Experimentally, it turns out that radiation forces from a cw-laser 
on a liquid-liquid interface are able to produce giant deformations (up 
to about 100 fJ,m), if the system is close to the critical point where 
the surface tension becomes small. We present a new model for such 
a fingerlike deformation, implying that the system is described as an 
optical fiber. One reason for introducing such a model is that the 
refractive index difference in modern experiments, such as those of the 
Bordeaux group, is small, of the same order as in practical fibers in 
optics. It is natural therefore, to adopt the hybrid HEn mode, known 
from fiber theory, as the fundamental mode for the liquid system. We 
show how the balance between hydrodynamical and radiation forces 
leads to a stable equilibrium point for the liquid column. Also, we 
calculate the narrowing of the column radius as the depth increases. 
Comparison with experimental results of the Bordeaux group yields 
quite satisfactory agreement as regards the column width. 

PACS numbers: 42.50.Wk, 47.61.-k, 42.25.Fx 

1 Introduction 

Force exerted by light on continuous matter has attracted attention since the 
days of Maxwell. In 1869 he was the first to predict and calculate the pressure 
from light reflecting off a surface; in 1901 the first attempts of experimental 
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verification were made by Lebedev, showing that the light pressure acts in- 
wards to a reflecting surface. Poynting pQ later extended the knowledge of 
radiation pressure by considering light incident from vacuum on the surface 
of a transparent dielectric, predicting the irradiation to produce an outward 
normal force irrespective of the angle of incidence. More recently, the ex- 
perimental paper of Ashkin and Dziedzic from 1973 is a pioneering work [2J. 
They illumined an air-water surface by a focused pulsed laser, verifying that 
the force was indeed outward directed. Another recent experiment along the 
same lines is that of Schroll et al. [3]. 

The radiation pressure problem is related to the well known Abraham- 
Minkowski controversy about the correct form of the electromagnetic energy- 
momentum tensor in a medium. That question has been discussed at vary- 
ing degrees of intensity ever since Abraham and Minkowski presented their 
energy-momentum expressions around 1910. In the simplest case, when the 
medium is isotropic and homogeneous, the difference turns up only in the ex- 
pressions for the momentum density g: in the Minkowski case g M = D x B, 
whereas in the Abraham case, g A =(l/c 2 )E x H, the latter expression satis- 
fying Planck's principle of inertia of energy g = S/c 2 , S being the Poynting 
vector. For an introduction to the Abraham- Minkowski problem, the reader 
may consult M0ller's book on relativity [I], or the review article [5] of one 
of the present authors. There is by now an extensive literature in this field; 
some papers are Refs. P El El El QUI EEH E321 E32]- Fortunately, for practi- 
cal purposes the difference between the Abraham and Minkowski predictions 
goes away in optics because the influence from the 'Abraham term' fluctuates 
out. The force can be calculated from the electromagnetic stress tensor parts 
only, and the stress tensors are equal in the two cases. We shall take this into 
account in the following, and simply call the force the Abraham-Minkowski 
(AM) force. 

A major reason why radiation forces has attracted increased interest in 
recent years is their practical usefulness. In biology and medicine the need 
of non-invasive techniques for manipulation of individual cells or complex 
molecules has led to the development of optical tweezers, invented by Ashkin 
et al. in 1986 [3]. Furthermore, fluid-interface instabilities driven by the 
relatively strong forces from electric field [T5l [T6] represent the cornerstone of 
many industrial processes such as electrospraying [17], ink-jet printing [IS] , 
and surface relief printing [T9] . 

An important progress in the application of pressure forces, is to make 
use of two liquids in the vicinity of the critical point. Then the surface ten- 
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sion is significantly diminished, and the effect of the pressure forces becomes 
enhanced. Traditionally, due to the competition with surface forces, the 
pressure effect has been rather minute. For instance, in the Ashkin-Dziedzic 
air-water pressure experiment [2j, with surface tension equal to a =73 mJ/m, 
the deflection of the water surface acted upon by a pulsed Nd YAG laser was 
only about 1 /mi. Working near the critical point, the deflection can be much 
higher, about 50 /j,m or more. Moreover, in a two-fluid system of surfactant- 
coated nano droplets in oil microemulsions the surface tension can be made 
more than one million times smaller than the usual air-water tension (cf, for 
instance, Ref. [20]). 

The main purpose of the present paper is to introduce a fiber model for the 
large finger-shaped deformations. Interesting work, both experimental and 
theoretical, has in recent years been done in this direction by the Bordeaux 
group, considering diverse effects arising from the focusing of laser beams 
onto a fluid-fluid interface having extremely low surface tension. Much of 
this work, up to 2002, is summarized in the PhD thesis of Casner [21]. There 
are several recent articles by this group [31 [221 [231 [2H [23 [22]- There is also 
a review paper [27] . and related papers of others such as [23]- Whereas the 
first experimental and theoretical investigations dealt with small deforma- 
tions in the hydrodynamic linear regime, the most recent experiments have 
demonstrated the appearance of giant surface deformations of order 50-100 
/im. Large finger-shaped structures have been observed, as well as liquid jets 
and shedding of microdroplets. 

A theoretical description of these finger-shaped deformations is still lack- 
ing. Nor will we in the present paper attempt to give a quantitative estimate 
of the large deformations. Instead, we will show one can calculate a stable 
equilibrium radius of the liquid column. As is conventional in fiber optics, 
we take the HE n mode to be dominant for a step-like fiber having a weak 
refractive index difference between core and cladding. This dominance oc- 
curs if the refractive index contrast is small, about 0.01. Actually, this fits 
quite well with the conditions of the Bordeaux experiments. It means that 
the central core becomes exposed to a laser intensity differing in shape and 
intensity from the conventional Gaussian intensity distribution. 

There is another advantage of this kind of wave description. The radia- 
tion pressure is predicted to be nonvanishing even on the vertical walls of the 
cavity. The magnitude of the pressure is of the same order as the pressure 
from a Gaussian beam on a flat surface at normal incidence. Such a pres- 
sure is not obtainable using a ray picture of the beam; in such the 
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(orthogonal) pressure on a parallel wall is simply zero. In our approach this 
problem is avoided. 

Section 2 gives, for the sake of readability, a brief account of the theoret- 
ical background in the linear case although this material is strictly speaking 
not new. In particular, we review how the equation of force equilibrium is 
solved when the input laser intensity is Gaussian. Section 3 then introduces 
our step-index fiber model. We use this model, as we indicated, to show how 
the balance of hydrodynamical and radiation forces leads to a stable equi- 
librium radius for a given power P in the laser beam. Comparison with the 
giant deformations observed in the Bordeaux experiments shows reasonably 
good agreement. 



2 Theoretical background, linear theory 

The electromagnetic force density f in an isotropic, nonconductive and non- 
magnetic medium (see, for instance, Refs. [3 [29]), is 



E 2 P 



' dn 2 
s dp 



tJ 



We employ SI units so that the relation e fiQ = 1/c 2 refers to a vacuum, and 
let e be a relative quantity so that the constitutive relations are D = e eE, 
B = /U H. The medium is assumed to be nondispersive. Only the first term 
in Eq. ([TJ) contributes in our case; as mentioned above this is the Abraham- 
Minkowski (AM) force 

f AM = _| j£; 2 Vn2 _ (2) 

The geometry is sketched in Fig. 1 (we assume illumination of the surface 
from below). The incident wave is taken to be monochromatic, = 
EW(r)e- <wt . The pi ane of incidence is formed by by the vectors kj and n; the 
angle of incidence is 9{ and the angle of transmission is 9t (see Fig. 1). Let 
E|| and Ej_ be the components of E parallel and perpendicular to the plane 
of incidence; then the respective transmission coefficients are 

sin 26 i sin 29 1 
11 = sin 2 (^ + ^)cos 2 (^-^)' 
sin 29i sin 29 t 

T± = sin^ + ^r (4) 
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Let I = eomciE® ) (averaged over the two polarizations) be the intensity of 
the incident beam. If a is the angle between and the plane of incidence, 
so cos a, sin a, we can write the surface pressure 

as [27J m 

U AM = _J_ n\-n\ cosOi 2 ^ + ^ ^ q + ^ ^ 
2c n 2 cosfc't 

Let now the liquid-liquid interface be represented asz = h(r, 9) in cylindrical 
coordinates (r = y/x 2 + y 2 ) and assume azimuthal symmetry, so that /i^ = 
= 0. Observing that cos^ = (1 + ^)" 1/2 , sin^ = h r (l + h 2 r )- 1/2 
{h r = dh/dr) as well as the corresponding relations for 6 t we obtain, upon 
substitution into Eq. 

2n x I\-N 

n = - 1 + N f( h r,<X), ( 6 ) 

where 

N = n 1 /n 2 < 1 
is the refractive index ratio and f(h r , a) is the function 

, (1 + N) 2 (. 2 1 + (3 - N 2 )h 2 + (2 - N 2 )ht 2 
f(h r , a) = , v ' < sin 2 a H . ' r y ' r cos 2 a 

[N+ y/l + h 2 {l - N 2 )] 2 { [Nh 2 + y/l + h^l-N 2 )} 2 

(7) 

Now it turns out that the polarization plays no important part in the present 
problem. This can be seen by plotting the normalized radiation force for 
different values of the angle a (not shown here; some more details can be 
found in Refs. [28] and [31]). Taking N = 0.986, a physically realistic value in 
the Bordeaux experiments, it is hard even to separate the curves for varying 
values for a. In the actual experiments the values of N were indeed near 
unity. This gives an expression for f(h r ,a) — > f(N,h r ), averaged over the 
angle a, that is more than accurate enough for our purpose: 



f(N h) = (1 + N ^ + (2-N 2 )h 2 + K + Nh 2 r S 
f[ ' r) l + ' (N + S) 2 (Nh 2 + S) 2 



(8) 



where S = ^Jl + h 2 (l - N 2 ). 

The system that we shall consider is comprised of two near-critical oil- 
emulsions separated by an interface, caused by the density contrast between 
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Figure 1: Laser wave incident upon an interface separating media 1 and 2, 
where the index of refraction ri\ < ri2- The electromagnetic surface force is 
directed towards the optically thinner medium irrespective of the laser light's 
direction of propagation. 

the liquids. It is useful to introduce the Bond-number, describing the relative 
strength of buoyancy in comparison to the Laplace force: 



Here ujq is the radius of the beam waist, and Iq is the capillary length, 



where a is the surface tension, p\ and p2 are the densities of the lower and 
upper liquids, and g is the gravitational acceleration. In the experiments of 
Casner and Delville, the Bond-numbers were in the range from 10~ 3 to about 
4. If B <C 1, the gravitational force is much weaker than the surface-tension 
force. 

An advantage of near-critical systems is the possibility to tune fluid prop- 
erties continuously by varying the temperature. Many physical quantities 
scale with temperature as oc [(T — Tc)/Tq]P, where f3 is a constant (see for 
instance [30]). Detailed information about liquid properties near to the crit- 
ical point is available elsewhere [23[ EU [27j EE]> an d will not be repeated 
here. 




(9) 




(10) 



Consider now the governing equation for the liquid-liquid surface eleva- 
tion. Under stationary conditions the elevation is determined by the balance 
of radiation pressure, surface tension, and gravity. Applying Laplace's for- 
mula for the pressure difference p 2 — Pi, we obtain [27] 12H] 



Apg ■ h(r) — 

r dr 



rh r 



l + hl 



2I(r) n 2 -n 1 

n x ■ f(N, h r 

c n 2 + n 1 



where Ap = p\ — p 2 > 0. The equation can be solved numerically once J(r) 
is known. 

Let us assume a Gaussian form: 



Kr) 



2P 



-2r 2 /w 2 



(12) 



where P is the total power. It is convenient to introduce the dimensionless 
variables R = r/uo, H(R) = h(r)/uo, whereby Eq. (ITT]) can be written in 
dimensionless form, 



BH--— 
RdR 



RH 



R 



-Fe 



-2R 2 



f(N,H R ) 



with 



2P\dn/dp\ 



(13) 



(14) 



ircguolc 

The function f(N, Hr) is the same as given by Eq. (jSJ) above, except from 
the replacement h r — > Hr. As before, B = (uo/lc) 2 is the Bond-number, 
and \dn/dp\ = 1.22 x 10~ 4 m 3 /kg j2TJ. The ratio N, and the capillary length 
la-, are the only temperature-dependent parameters. Thus, there are three 
parameters in Eq. (|13[) . namely the temperature difference AT, the beam 
waist uq, and the beam power P. 

Numerical solutions of Eq. (jl~3l) are shown in Fig. [2j Similar results are 
presented in Ref. [28], except that in the present case the theoretical curve 
is compared directly with experimental data from the Bordeaux group. Pro- 
files of the deformations are extracted from images taken by a CCD-camera 
(provided by Jean-Pierre Delville, personal communication). As is seen, the 
numerical solution gives excellent results for a laser power of P=450 mW, 
beam waist ujq = 5.3 /im, and AT = T — Tc = 3.5 K. If the laser power 
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Theory, P = 450 mW 

* Experiment, P = 450 mW 
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Figure 2: Experimental versus numerical results at temperature such th 
AT = T — Tq = 3.5 K, and with a laser beam waist of loq = 5.3 /xm at 
power of 450 mW. 
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is increased, the correspondence between theory and experiment becomes 
poorer. For instance, insertion of the large power P = 1200 mW in the for- 
mulas would lead to a considerable overprediction (about 50%) of the surface 
displacement. 

This concludes our overview of the linear theory. 

3 Electromagnetic wave analysis 

When the power of the laser is increased (for a fixed beam width) beyond a 
certain threshold, the behavior of the liquid-liquid interface deformation be- 
comes nonlinear. It has been observed that a "shoulder" appears in the case 
where light propagates in the upward direction (cf., for instance, Ref. [2~Tj). 
If the light propagates downwards, the deformations become very large and 
cylinderlike, eventually forming a jet. The deformations may even form a 
liquid bridge crossing a central layer of fluid. In such it is natural 

to suggest that the deformations, or rather the liquid column, may be re- 
garded as an optical fiber, capable of guiding the laser-light due to internal 
reflections. As already mentioned, this is the main new idea of the present 
paper. Furthermore, the deformations and columns tend to have vertical 
walls, which corresponds to \h r \ — > oo. 

Examining the expression ([2D for the radiation pressure, it can be seen 
that (see Fig. 3) 

n AM -> 0, when |h r | -> oo, 

implying that the radiation force vanishes, leaving only the forces of surface- 
tension and gravity, which both act against any deformation of the interface. 
This result is obviously unphysical, since both steep- walled deformations, and 
the liquid columns, are stationary structures. This means that the following 
two usual assumptions: 

(i) The laser-beam intensity distribution is Gaussian, and stays Gaussian 
also in the region being deformed; 

(ii) ray optics may be used to describe the path of a light beam within 
the deformation, 

must be invalid. Regarding the intensity distribution, one can argue that the 
geometrical change of the illuminated interface will result in new boundary 
conditions that the electromagnetic field must satisfy. In turn, this gives rise 
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Figure 3: The function f(N, h r ) goes to zero as h r grows large, i.e. when the 
walls of the deformation, or the liquid column, are (nearly) vertical. The only 
remaining forces is then the Laplace force and the hydrostatic pressure-force, 
both acting inwards. This implies that the column or deformation should 
collapse in on itself, which it evidently does not. 

to a different, and unknown, intensity distribution within the huge defor- 
mations characteristic of the non-linear regime. Also, ray-optics may give 
a poor description of the propagation of the electromagnetic waves within 
the structure. In our optical fiber model the core of radius r = a is the 
upper liquid, and the cladding is the surrounding liquid (the lower liquid). 
Solving Maxwell's equations for such a geometry, and relating the power flow 
through the structure to the power of the incident laser beam, one should be 
able to obtain a more reliable expression for the radiation pressure, and give 
an estimate for its magnitude. 

3.1 Derivation of the modified radiation force 

Assume as before plane wave propagation and sinusoidal time-dependence of 
the fields, E, B oc e~ wt . Maxwell's equations may be combined to give the 
following Helmholtz wave equations for H and E (cf., for instance, Ref. |33j): 

V 2 H + ^fH = M„(Vn 2 )xE, (15) 
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Figure 4: Section of the deformation modelled as an optical step-index fiber. 
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(16) 



By taking the z-components of Eqs. (I15p and ({TBI) , and eliminating the trans- 
verse field components, one finds the two-dimensional scalar wave equations 
for the system, 



CO 



V 2 t H z + 1 2 H z - — (V t n 2 )-V t H : 



and 



7c 



V/E z + YH z -\— (V t n')-V t E 



■^■[V^xVA], (17) 

T 



^■[V^xVA], (18) 



7*n 



where V 2 = V 2 — d 2 is the transverse part of the Laplacian operator. The 
solutions to these equations are 



{E z ,H z } = {A e} A h }J m (iry 



imO+ikzZ—iiiit 



and 



{E z ,H z } = {B e ,B h }K m ((3r)e" 



imd- s rih z z—iuit 



r < a 



r > a. 



(19) 
(20) 
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Here, m is an integer, J m is the m'th Bessel function of the first kind, and 
K m is a modified Bessel function. The radius of the deformation (the core 
radius) is r = a. The constants A e , Ah, B e , and Bh are determined from the 
boundary conditions at r = a. Furthermore, 

7 2 = co 2 n 2 2 /c 2 - kl = klnl - kl = k\ - kl (r < a) 



(3 2 = k\ - u 2 n\/c 2 = k 2 z - kln\ = k 2 z - kf, (r > a) 

are the radial propagation constants. k z is the longitudinal wave number, 
given by k z = Accost?, with •& being the angle of propagation. See Fig. 0J 
Note that ko = uo/c refers to the vacuum, k2 = ^o n 2 refers to the core, and 
ki = korii refers to the cladding. 

Inserting the solutions (|T9|) and (|20|) into the wave equations (IT7|) and 
( TT8|) . one can find the radial and azimuthal field-components. Detailed cal- 
culations are given in the books of Stratton [2H] arid Okamoto [31] • One has 
to observe that there are multiple solutions to Maxwell's equations for the 
step-index geometry, corresponding to different modes of propagation. We 
restrict ourselves here to giving the expressions for the radial and azimuthal 
components for the electric field inside the fiber, r < a, 

E r = U^J' m M - ^^J m ( 7 r)) 
I 7 TT J 

E 9 = - { k -^J m M + i^^J'^r)) e^—', (21) 

E z = A e J m Me me+lk * z -^, 
and similarly the if -field for r < a, 

I 7 l 2 r J 

He = l-^^JrnM+i^^JLM) (22) 
I 7 2 r 7 J 



H z = A h J m (jr)e 



imd+ik z z—iujt 
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The prime on the Bessel functions indicates differentiation with respect to 
the argument jr. 

The aim now is to determine the allowed discrete angles at which the 
light rays may propagate, corresponding to the allowed modes of propagation. 
We begin by defining the normalized transverse wave numbers as 



u = 7a 



aJB - k 2 



(23) 



w = (5a = ay/kl - k\. (24) 
Furthermore, the wave numbers u and w are related, from Eq. fl36|) . as 

u 2 + w 2 = k\{n\ - n\)a 2 = v 2 , (25) 

referred to as the normalized frequency. The boundary conditions are that 
the tangential field components are continuous across r = a. From them the 
general dispersion relation is constructed, valid for all values of N — ni/n 2 < 
1 



J' m {u) + K'Jw) 



uJm.(u) wK m (w) 



J ' m ( U ) +N 2 K 'm( W ) 



Uj m (u) 



wK m (w) 



m 



1 1 

u 2 w 2 



1_ 



w 



(26) 



3.2 Fundamental mode of the step-index fiber. Calcu- 
lation of the equilibrium radius 

To simplify the discussion, we shall from now on consider only the funda- 
mental mode of the step-index optical fiber. This mode corresponds to the 
m = 1 solution of the scalar wave equation, and is called a hybrid mode. 
The fundamental mode of a step-index optical fiber is the HE m=1 i=1 mode, 
corresponding to both E z and H z nonzero [3U ES] , in contrast to the trans- 
verse electric (TE) and transverse magnetic (TM) modes which correspond 
tom = having E z and H z equal to zero respectively, but which possess low 
cutoff frequencies. 

The index I = 1 corresponds to the first root of the dispersion relation 
satisfying k\ < k z < k 2 . The HE n mode has no cutoff frequency; this 
is the reason why it is regarded as the fundamental mode of a step-index 
fiber. For a typical fiber in optics, the relative refractive index difference 
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A = 1 — n\jni is of the order of 0.01. This suits well with the conditions of 
the Bordeaux experiments: with the typical value of N = 0.986 as mentioned 
above, we get A = 0.014, which is quite appropriate. Under such conditions 
the confinement of light in the core is not so tight (this is called the weakly 
guiding approximation). The approximation N ~ 1 allows one to simplify 
the theory of optical fibers very much, and to obtain clear-cut results. The 
accuracy of adopting the HE n mode as a basis in our model should be more 
than sufficient. 

A sketch of the intensity distribution for this mode is given in Fig. 5. 
Numerically solving the dispersion relation (|26|) for a given fiber radius a, 
we obtain the wave number k z corresponding to that particular radius. Re- 
peating this procedure for a range of radii, and with all field components 
being known, we can calculate the components of the Abraham-Minkowski 
surface force density on the interface of a fiber with a given diameter. The 
expression for the AM surface pressure reads, 



El + El + N~ 2 E 2 r 



n. 



Here 



n=(l + h 2 x + h 2 y y 1 / 2 {h x ,h y ,-l) 



(27) 



(28) 



is the normal vector to the interface, pointing from the internal medium 2 to 
the external medium 1. If S = |(E x H*) denotes the Poynting vector, the 
total power P carried by the optical fiber plus the cladding is given by 



P= / S z rdrd6 = - / (E r H* e - E e H*)r dr d9 . 
Jo Jo 2 Jo Jo 

Some calculation yields 



1 r 2n 



TV 



P — Pcore + Pclad ~ ^0 ca2 \^\ 2 



nzFi^Jj+n^GiKo,^) 



u 



where 



F( J , Ji) = 



and 



G{K Q ,K X ) = 



k 2 n 2 r 9 



'l + s 2 )(Jl + Jl)--(l + sfJl 



1 1 



(29) 



(30) 



(31) 



^\(l + s 2 )(K!-K 2 ) + ^(l + s) 2 K! 
w A L w z 



+ -K K 1 . (32) 
w 
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Here, it is understood that J n , K n = J n (u), K n (u). The constant A is equal to 
A e above, and B e has been eliminated via the boundary condition E z (a—) = 
E z (a+). The parameter s is given by 



\M Z W J 



uj\(u) wKi(w) 



(33) 



(cf. Ref. [M]). When the steepness is large (h r ^> 1) the expression for the 
Laplace force in the radial direction takes on the simple form 

h = ~ lim °*^L= = (34) 



The force balance may then be written as 

- Apgh(r = a) + - = n AM (r = a), (35) 



This equation can be interpreted physically: The right hand side is the AM 
surface force density acting outwards. It is positive because n 2 > n\. This 
force is balanced by the Laplace force a/a acting inwards, plus the net hy- 
drostatic pressure —Apgh(r = a) = — (pi — p 2 )gh(r = a) which also acts 
inwards because pi > p 2 . Notice that h(r) is a negative quantity. We do 
not regard a as a fixed parameter, but rather determine the magnitude of 
the radiation pressure for a range of given values of a, and subsequently 
find the corresponding equilibrium radii where the radiation pressure exactly 
balances the forces of surface tension and hydrostatic pressure. Ergo, the 
thickness of the liquid column will vary with depth. For a given depth, the 
hydrostatic pressure difference is known, and then it is easy to determine 
the new equilibrium radius a. Notice that we still regard the structure as an 
ideal step- index fiber, although in fact it is not. The fiber rather resembles a 
tapered cylinder, which would require a different solution to Maxwell's equa- 
tions to get exact results. But in our model, it seems reasonable to regard 
the fiber locally (with depth) as an optical cylindrical guide with zero taper. 



Figure 6 shows how the electromagnetic radiation AM force varies versus 
the fiber radius a , for given values of AT and P. The combined pressure of 
surface tension (the Laplace force) and gravity (hydrostatic pressure) is also 
plotted. Intersect points correspond to equilibrium radii. Notice the middle 
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Figure 5: The intensity distribution for the HEn-mode in a step-index fiber 
of radius a = 2 /im, at a laser beam power of 500 mW, and with AT = 3.5 K. 
Broken lines indicate the fiber boundary. 
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set of curves with power P = 150 mW. Here, the termination point of the 
liquid column can clearly be seen. As the depth increases beyond 550 fim, 
the outward acting AM force can no longer support the liquid column due 
to increasing hydrostatic pressure and surface tension. The model therefore 
predicts that at these parameter values, the liquid column can maximally 
reach a depth of ~ 500 /jjd. before collapsing, in good correspondence with 
observed break-ups of the liquid columns into droplet spraying jets at depths 
of around 1 jum. In our simple model we do not take into account the fluid 
flow inside the column; nevertheless, the model offers a simple and intuitive 
explanation of the observed pinch off of the columns. 

With increasing temperature difference AT the radiation force becomes 
weaker relative to the Laplace force and the hydrostatic pressure (details not 
shown here). This behavior is expected, since the surface tension is tempera- 
ture dependent and an increase in AT leads to a higher surface tension. Now 
the hydrostatic pressure slightly increases as the density contrast increases, 
but this is not of great importance since the term Apgh is very small. The 
net behavior is that for higher AT the equilibrium radii become smaller. 

In Fig. 7, the two left panels refer to temperature differences, AT = 
{3.5, 1.5} K, and to a moderate power, P =200 mW. These parameter values 
make our model quite justifiable, as can be seen from the following argument: 

• When AT is small the relative refractive index becomes also small, in 
our case amounting to 0.55% (i.e., A = 1 — ni/n 2 = 0.0055). The 
microemulsion liquid used in the Bordeaux experiments had a temper- 
ature dependent refractive index of around 1.46 ([21]). The laser used 
was an argon laser operating in the TEMoo mode, with a vacuum wave- 
length of Ao = 0.514 /im. From these values we calculate the normalized 
frequency 



• The second point is that for moderate powers the equilibrium radii for 
the liquid column become small. Physically, the reason is that the AM 
outward pressure is not strong enough to widen the column very much. 
The two left panels in Fig. 7 show that for P = 200 mW the radii are 
about 1 /im. It is instructive here to compare with Fig. 3 in Okamoto's 
book [34J, which shows the occurrences of multiple modes in a typical 
single-step fiber. As long as v is less than about 2.4, only the basic HE n 
mode appears. For increasing v, up to 3.5, perhaps even higher, there 




(36) 
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occurs a mixture with the next HE 21 mode, but we find it reasonable 
to expect that the basic mode is the dominant one at least up to about 
v = 3. Insertion of a = 1 /im in Eq. fl36l) yields v ~ 2. Accordingly, 
the single-mode model ought in this case to be a viable one. Also, we 
note from Fig. 7 that the uniform- cylinder model of the column turns 
out to be a good approximation, as the slope of the column is large, 
\h r \ ~ 10 4 . 

Actually, it seems that the above limits for the applicability of the single- 
mode model may be stretched considerably. The third panel of Fig. 7 shows 
the predicted column width when the power is high, P = 1750 mW. This is 
in violation of the single-mode condition, but still the calculated results yield 
a figure visually quite similar to the experimental deformation shown in the 
rightmost frame. The last-mentioned frame, taken from Ref. [2T], shows a 
long jet with large radius, a « 7 /im. The good agreement between theory 
and experiment indicates that the single-step model is quite robust. 

The last point that we shall focus attention on, is the observed termina- 
tion, or "pinch off", of the liquid column at large depths where the radius 
becomes small. We can obtain a semi-quantitative explanation of this effect 
on the basis of our model. Namely, when a becomes smaller than about 
0.6 /im the repulsive AM radiation force is no longer able to withstand the 
inward directed Laplace and hydrostatic forces (there is no equilibrium point 
in Fig. 6), and the column has to collapse. What one actually observes in 
experiments, is that the column breaks up into a row of small droplets. In 
Fig. 7 (a) and (b), the theoretical pinch-off of the columns is seen to occur 
at depths of around 1000 /im. At equal power, the column stretches deeper 
when the surface tension is lower (A T =1.5 K). It is moreover observed that 
the break-up occurs in the case of illumination from above only. Our model 
is however unable to explain that particular effect. No theory for it seems to 
be known. 
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Figure 6: Color online. The electromagnetic radiation pressure (AM force) 
plotted as a function of varying fiber radius a, at given temperature AT and 
beam power P. 
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4 Summary 



Let us summarize our work as follows: 

(1) Our main idea has been to model the large fingerlike deformations, 
the liquid columns and liquid jets seen in laser optics experiments, as optical 
waveguides or fibers. We have identified the electric field components of a 
step-index fiber, identifying the HE n mode as the fundamental mode. The 
HE n mode is a hybrid mode, thus more complicated than the conventional 
TE and TM modes in that the axial electromagnetic fields E z and H z are not 
zero, but it is precisely this mode that is the dominant one in conventional 
step-index fibers when the index contrast iV = n\jni is close to 1 [31]. The 
fact that the giant deformations are seen experimentally in the vicinity of 
the critical point when N ~ 1, makes our fiber model very natural from a 
physical viewpoint. Moreover, the experiments performed by the Bordeaux 
group [2H [31 [221 [231 [211 [251 [26] correspond just to N « 1. 

(2) In Sect. 2 we reviewed the linear theory, showing how the elevation of 
the interface can be calculated with quite good accuracy for low laser powers 
P when comparing with the Bordeaux experiments. No attempt was how- 
ever made to calculate the elevations in the nonlinear, high-power case. In 
particular, we did not consider the formation of a "shoulder" as seen exper- 
imentally upon illumination from below. No theory for the nonlinear case 
seems so far to exist. Instead, we made use of our model to show how the 
balance between hydrodynamical and radiation forces leads to the establish- 
ment of stable radii for the liquid column. See Fig. 6. This means that a 
physical explanation is given for the establishment of these giant structures. 
We emphasize that such a result is not obtainable from conventional radia- 
tion theory assuming a Gaussian intensity profile for the incident laser beam. 
Our intensity distribution described by the HE n mode is different from the 
Gaussian form. See the more detailed discussion on this point in the first 
part of Sect. 3. 

(3) Our comparisons with experiments are all made with reference to the 
Bordeaux group. Figure 7 shows how our calculated values for the column 
radius a correspond quite well with the values observed. In particular, the 
narrowing of the width for increasing depths is reproduced. Such a narrowing, 
by a factor of 2 or more, is typically seen in experiments (Jean-Pierre Delville, 
personal communication). 
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Figure 7: Color online. Calculated interface deformation resulting from the 
balance between the modified AM-force, the Laplace force, and the hydro- 
static force. Panel (a) and (b) have P = 200, and AT = {3.5, 1.5}K, re- 
spectively. At this moderate power level and values of AT, the radii of the 
columns are, according to the model, around 1 ^m, thus satisfying the single- 
mode condition of the step-index fiber. At high power P = 1750 mW and 
AT = 6K, shown in panel (c), the column is predicted to be slightly wider 
with a radius from 1.5 to 3 /im. The picture (d) is taken from Ref. [2T], and 
shows the liquid column resulting from illuminating the surface from above 
with a laser at P = 1750 mW, and with T — Tc = 6 K. In Ref. [21] it is stated 
that the jet is about 1000 /im deep, with a radius around 7/xm, i.e., of the 
same order of magnitude as that calculated in our model. 
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Note added 

After the above work was completed, we have become aware of two recent 
papers of the Bordeaux group [361 [37] . These papers report on both theo- 
retical and experimental work. As for the theoretical parts, they are closely 
related to our approach above. 
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